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Abstract 

The existence of a Noether symmetry for a given minisuperspace cosmolog- 
ical model is a sort of selection rule to recover classical behaviours in cosmic 
evolution since oscillatory regimes for the wave function of the universe come 
out. The so called Hartle criterion to select correlated regions in the configu- 
ration space of dynamical variables can be directly connected to the presence 
of a Noether symmetry and we show that such a statement works for generic 
extended theories of gravity in the framework of minisuperspace approxima- 
tion. Examples and exact cosmological solutions are given for nonminimally 
coupled and higher-order theories. 
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I. INTRODUCTION 



Several points of view can be adopted in order to define quantum cosmology. It can 
be considered as the first step toward the construction of a full theory of quantum gravity. 
Besides, it concerns finding initial conditions from which our classical universe is started. 
However, with respect to other theories of physics as electromagnetism, general relativity or 
ordinary quantum mechanics, boundary conditions for the evolution of the system "universe" 
cannot be set from outside. There we need a fundamental dynamical law (e.g. Maxwell's 
or Einstein's equations or Schrodinger's equation) and then we impose, from the outside, 
the initial conditions. In cosmology, by definition, there is no rest of the universe so that 
boundary conditions must be a fundamental law of physics. In this sense, a part the fact that 
quantum cosmology is a workable scheme to achieve quantum gravity, it can be considered 
as an autonomous branch of physics due to the issue of finding initial conditions |L|. 

However, not only the conceptual difficulties, but also mathematical ones make quantum 
cosmology hard to handle. For example, the superspace of geometrodynamics || has infinite 
degrees of freedom so that it is practically impossible to integrate the full Wheeler-DeWitt 
(WDW) equation. Furthermore, a Hilbert space of states describing the universe is not 
available ||] . Finally, it is not clear how to interpret the solutions of WDW equation in the 
framework of probability theory. Several interpretative schemes have been proposed but the 
concepts of probability and unitarity are in any case approximate. Their validity is limited 
by the accuracy of the semiclassical approximation and strictly depends on the suitable 
definition of probability current ||, [|J. 

Despite these still unsolved shortcomings, several positive results have been obtained and 
quantum cosmology has become a sort of paradigm in theoretical physics researches. For ex- 
ample the infinite-dimensional superspace can be restricted to opportune finite-dimensional 
configuration spaces called minisuper spaces. In this case, the above mathematical difficul- 
ties can be avoided and the WDW equation can be integrated. The so called no boundary 
condition by Harte and Hawking || and the tunneling from nothing by Vilenkin Q give 
reasonable laws for initial conditions from which our classical universe could be started. 

The Hartle criterion || is an interpretative scheme for the solutions of the WDW equa- 
tion. Hartle proposed to look for peaks of the wave function of the universe: If it is strongly 
peaked, we have correlations among the geometrical and matter degrees of freedom; if it is 
not peaked, correlations are lost. In the first case, the emergence of classical trajectories 
(i.e. universes) is expected. The analogy to the non-relativistic quantum mechanics is im- 
mediate. If we have a potential barrier and a wave function, solution of the Schrodinger 
equation, we have an oscillatory regime on and outside the barrier; we have a decreasing 
exponential behaviour under the barrier. The system behaves classically in the oscillatory 
regime while it does not in the exponential case. The situation is analogous in quantum 
cosmology: Now the potential barrier has to be replaced by the superpotential U(h l \(p), 
where are the components of the three-metric of geometrodynamics and (p is a generic 
scalar field describing the matter content. More precisely, the wave function of the universe 
can be written as 



where m p is the Planck mass and S is an action. A state with classical correlations must be 
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a superposition of states of the form This type of state can be expressed as a coherent 

superposition of eigenstates of operators that commute with the constraints and correspond 
to constants of the motion. A superposition of this kind can be approximated by a WKB 
state where 

S = S + m P 2 S 1 + 0(mp 4 ) , (1.2) 

is the expansion of the action. We have to note that there is no normalization factor due 
to the lack of a probability interpretative scheme. However, in this approximation, it is 
possible to define a localized prefactor as shown in and we can define a quasiclassical 
state to describe an approximate classical behaviour and a semiclassical one to describe a 
product of a part that is quasiclassical and a part that is not. This quasiclassical state can 
be a coherent one if it is a superposition of states in the sense discussed in 0. Transition 
amplitudes of istantaneous eigenstates are discussed in |8j . 

Considering the action (|1.2j ) and inserting it into the WDW equation and equating similar 
power of m p , one obtains the Hamilton- Jacobi equation for Sq. Similarly, one gets equations 
for Si, S 2 , ■ ■ ., which can be solved considering results of previous orders. We need only S 
to recover the semi-classical limit of quantum cosmology |J. If So is a real number, we get 
oscillating WKB modes and \I/ is peaked on a phase-space region defined by 

where 7Ty and tt^ are classical momenta conjugates to and tp. The semi-classical region of 
superspace, where \& has an oscillating structure, is the Lorentz one otherwise it is Euclidean. 
In the latter case, we have S = il and 

y~ e - m i I , (1.4) 

where I is the action for the Euclidean solutions of classical field equations (istantons) . Given 
an action 5*0, Eqs. (|1.3j ) imply n free parameters (one for each dimension of the configuration 
space Q = {h %3 , <p}) and then n first integrals of motion. However the general solution of the 
field equations involves 2n — 1 parameters (one for each Hamilton equation of motion except 
the energy constraint). Consequently, the wave function oscillates on a subset of the general 
solution. In this sense, the boundary conditions on the wave function (e.g. Harte-Hawking, 
Vilenkin, others) imply initial conditions for the classical solutions. 

To be more precise, if the wave function \1/ is sufficiently peaked about some region in 
configuration space, we predict that we will observe the correlations between observables 
which characterize this region. If \1/ is small in some region, we predict that observations of 
the correlations which characterize this region are precluded. Where \I> is neither small nor 
sufficiently peaked, we do not predict anything. 

For example, given the measured value of the Hubble constant and mass density, we 
would like a "good" wave function for the universe to be peaked around a distribution of 
galaxies consistent with that which is observed. It is crucial to recognize that the wave 
function does not predict a specific value for H , or specific locations for the galaxies, but 
rather a "correlation" between these observables. 
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Halliwell has shown that an oscillatory wave function of the form ( |1 . 1| ) predicts a 
correlation between the canonical coordinate q and the momentum Ti q of the above form ix q = 
m 2 P dS/dq. In other words, taking into account minisuperspace models, if oscillatory regimes 
of the WDW wave function exist, we are able to recover correlations among variables and 
then classical behaviours emerge. Again, in the classically allowed region, the semiclassical 
approximation to the WDW wave function yields just such oscillatory solutions. 

A simple minisuperspace example can be constructed. The ansatz for the wave function 

is 



D im'pS(a) 



;i.5) 



where the canonical variables q coincides to the scale factor of the universe a and the phase 
is a slowly varying function of the scale factor. Since we wish to investigate the classical 
limit m 2 p — > oo (which corresponds to h —>■ of ordinary quantum mechanics), we use the 
expansion (|1.2|) . Inserting such an ansatz into the WDW equation 



deduced by the action 



da 2 \2GJ 



V(a) = 



;i.6) 



A 



3tt 
AG 



dt 



-a a + a 1 



;i-7) 



of a Friedman-Robert son- Walker (FRW) closed universe (a is a constant) with the canon- 
ical momentum ir a given by 



7tn 



3tt 
2G l 



-aa . 



;i.s) 



we get a set of differential equations, one for any order of rrip , which as we said above, 
can be solved sequentially. The semiclassical approximation to the wave function obtains 
by working only to first order. We get 



Sn 



da' t 



/3tt\ 2 2 /a 4 a 2N 



;i.9) 



and 



Si 



2 I da 



;i.io) 



Thus the oscillatory semiclassical wave function \^ oc exp(z5*o) is peaked about a region 
of minisuperspace (every point of which represents a closed FRW model) in which the 
correlation between the coordinate and momentum (scale factor and expansion rate), rr a = 
dS /da, holds good. 



Using Eq. ( |1.8|) for n a , the correlation reduces to 
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la 2 



a -l, (1.11) 




which is nothing else but the (0, 0) Einstein equation for a FRW spacetime. If 

a = ^|, A = 8ttG/w, (1.12) 

where p vac is a constant density, we get the solution 

a(t) = ao cosh(ciQ 1 t) , (1.13) 

which is an inflationary behaviour for a closed FRW model. Thus, by this simple example, in 
the region of minisuperspace where the wave function oscillates, a classical FRW spacetime, 
obeying the (classical) Einstein equation emerges. 

The issue is now if there exists some method capable of selecting such constants of motion 
which, being first integrals of motion, allow to find correlations between classical variables 
and conjugate momenta in minisuperspace models. In other words, can the emergence of 
classical trajectories be implemented by some general approach without arbitrarily choosing 
regions of the phase-space where momenta ( |1.3| ) are constant? Achieving this result means to 
obtain oscillatory subsets of WDW wave function where one gets correlations. Consequently 
classical regime are recovered and the Hartle criterion holds (at least in the framework of 
the minisuperspace approximation). For the full theory, i.e. without considering simple 
minisuperspace models, the issue is more delicate since we have to ask for superpositions of 



the form ( |1 . 1| ) which yield peaked wave packets so that the Hartle criterion holds 0. In this 
case, also the issue of universe "creation", as a particle creation problem, has to be faced 
considering the way in which the quantum-classical transition is achieved ||. 

In this paper, we want to restrict to a more specific (and simple) question. We want to 
show, for general extended gravity minisuperspace models, that the existence of a Noether 
symmetry implies, at least, a subset of the general solution of the WDW equation where 
the oscillating behaviour is recovered. Viceversa, the presence of a Noether symmetry gives 
rise to the emergence of classical trajectories. This analysis is performed in the context 
of the minisuperspace approximation and, for classical trajectories, we mean solutions of 
the ordinary Einstein equations. The existence of a Noether symmetry for a dynamical 
model is a general criterion to search for constants (first integrals) of motion so that, given a 
minisuperspace model exhibiting such symmetries we obtain certainly correlations and then 
classical behaviours. This statement, in our knowledge, has never been done also if constants 
of motion have been systematically used in quantum cosmology since at least fifteen years. 

The layout of the paper is the following. Sect. 2 is devoted to the Noether Symmetry 
Approach and to its connection to quantum cosmology. In Sect. 3, we apply the method to 
minisuperspace models of nonminimally coupled theories of gravity, while the same is done 
in Sects. 4 and 5 for higher-order theories. Discussion and conclusions are drawn in Sect. 
6. 



5 



II. THE NOETHER SYMMETRY APPROACH AND QUANTUM COSMOLOGY 



Minisuperspaces are restrictions of the superspace of geometrodynamics. They are finite- 
dimensional configuration spaces on which point-like Lagrangians can be defined. Cosmo- 
logical models of physical interest can be defined on such minisuperspaces (e.g. Bianchi 
models). 

Before taking into account specific models, let us remind some properties of the Lie 
derivative and the derivation of the Noether theorem [|K]]. Let L X be the Lie derivative 

{L x u)i = juigiO , (2.1) 

where uj is a differential form of lZ n defined on the vector field £, g\ is the differential of 
the phase flux {g t } given by the vector field X on a differential manifold M.. Let p t = p g ~t 
be the action of a one-parameter group able to act on functions, vectors and forms on the 
vector spaces C°°(M), D(A4), and A(A4) constructed starting from M.. If g t takes the 
point m G M in gt(m), then p t takes back on m the vectors and the forms defined on gt{m); 
p t is a pull back [jO]]. Then the property 

p t+s = ptp s (2.2) 

holds since 

9t+ s = g t °g s - (2.3) 
On the functions /, g G C°°(J\A) we have 

Pt(fg) = (p t f)(Ptg); (2.4) 

on the vectors X, Y G V(M), 

p t [X,Y] = [p t X,p t Y]; (2.5) 

on the forms w,|iG A(Ai) 

p t {u A p) = (p t u) A (p t n) . (2.6) 

Lx is the infinitesimal generator of the one-parameter group p t , and, being a derivative on 
the algebras C°°(Ai), V(Ai), and A(A4), the following properties have to hold 

L x (fg) = (L x f)g + f(L x g), (2.7) 
L x % Z] = [L X Y, Z] + % L X Z] , (2.8) 
L x (ujAfi) = (L x u) Afi + u A(L X fi), (2.9) 

which are nothing else but the Leibniz rules for functions, vectors and differential forms, 
respectively. Furthermore, 
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Lxf = Xf, (2.10) 
L X Y = adX{Y) = [X, Y] , (2.11) 
L x du = dL x u , (2.12) 

where ad is the self-adjoint operator and d is the external derivative by which a p-form 
becomes a (p + l)-form. 

The discussion can be specified by considering a Lagrangian £ which is a function defined 
on the tangent space of configurations TQ = <ji}. In this case, the vector field X is 

where dot means derivative with respect to t, and 

L x C = XC = a\ q ) d ^ i+ a\ q ) d ^. (2.14) 

The condition 

L x C = (2.15) 

implies that the phase flux is conserved along X: This means that a constant of motion 
exists for C and the Noether theorem holds. In fact, taking into account the Euler-Lagrange 
equations 

d dC dC i * 

TtW-W = °' (2 ' 16) 



it is easy to show that 



If (P15D holds, 



(217) 



20 = ^ (2.18) 



is a constant of motion. Alternatively, using the Cartan one-form 

dC 

Oc = ^*t (2.19) 

and defining the inner derivative 

ix0c=<Oc,X>, (2.20) 

we get, as above, 

ixBc = So (2.21) 
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if condition ( 2.15Q holds. This representation is useful to identify cyclic variables. Using a 
point transformation on vector field ( [2.13D , it is possible to getQ 



X = {i x dQ k 







dQ k 



+ 



d_ 

dt 



ixdQ h 







(2.22) 



If X is a symmetry also X has this property, then it is always possible to choose a coordinate 
transformation so that 



i x dQ l = 1 , i x dQ i = , i ^ 1 



and then 



X 



d 



dC 



W dQ 1 



0. 



(2.23) 



(2.24) 



It is evident that Q l is the cyclic coordinate and the dynamics can be reduced However, 
the change of coordinates is not unique and a clever choice is always important. Furthermore, 
it is possible that more symmetries are found. In this case more cyclic variables exists. For 
example, if Xi,X 2 are the Noether vector fields and they commute, [Jfi,X 2 ] = 0, we obtain 
two cyclic coordinates by solving the system 



(2.25) 



i Xl dQ i = Q, Ml; ix 2 dQ i = 0, i^2. 



If they do not commute, this procedure does not work since commutation relations are 
preserved by diffeomorphisms. In this case 



X3 — [Xi , X 2 



is again a symmetry since 



(2.26) 



(2.27) 



Lx 3 £ — Lx ± Lx 2 £ — Lx 2 Lx ± £ — . 

If X 3 is independent of X\ , X 2 we can go on until the vector fields close the Lie algebra 

A reduction procedure by cyclic coordinates can be implemented in three steps: i) we 
choose a symmetry and obtain new coordinates as above. After this first reduction, we 
get a new Lagrangian L with a cyclic coordinate; ii) we search for new symmetries in this 
new space and apply the reduction technique until it is possible; iii) the process stops if 
we select a pure kinetic Lagrangian where all coordinates are cyclic. This case is not very 



1 We shall indicate the quantities as Lagrangians and vector fields with a tilde if the non-degenerate 
transformation 

Q l = Q\q), Q%) = ^f 

is performed. However the Jacobian determinant J =|| dQ % /dq J \\ has to be non-zero. 
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common and often it is not physically relevant. Going back to the point of view interesting 
in quantum cosmology, any symmetry selects a constant conjugate momentum since, by the 
Euler-Lagrange equations 

^ = 0^^ = S,. (2.28) 
dQ* dQi y 1 

Viceversa, the existence of a constant conjugate momentum means that a cyclic variable has 
to exist. In other words, a Noether symmetry exists. 

Further remarks on the form of the Lagrangian £ are necessary at this point. We shall 
take into account time-independent, non-degenerate Lagrangians £ = i.e. 

d£ d 2 £ 

- = 0, det/f^detll^ll/O, (2.29) 

where is the Hessian. As in usual analytic mechanics, £ can be set in the form 

£ = T(q\?)- V(q<), (2.30) 

where T is a positive-defined quadratic form in the qi and V(q l ) is a potential term. The 
energy function associated with £ is 

d£ 

E c = —c?-£{q^q 3 ) (2.31) 
and by the Legendre transformations 

d£ 

H = Ktf-£{q 3 ,qi), ^ = — , (2.32) 

we get the Hamiltonian function and the conjugate momenta. 

Considering again the symmetry, the condition ( [2.15| ) and the vector field X in Eq.( j2.13 ) 
give a homogeneous polynomial of second degree in the velocities plus an inhomogeneous 
term in the qK Due to ( |2.15| ), such a polynomial has to be identically zero and then each 
coefficient must be independently zero. If n is the dimension of the configuration space 
(i.e. the dimension of the minisuperspace) , we get {1 + n(n + l)/2} partial differential 
equations whose solutions assign the symmetry, as we shall see below. Such a symmetry is 
over-determined and, if a solution exists, it is expressed in terms of integration constants 
instead of boundary conditions. 

In the Hamiltonian formalism, we have 

[E j ,H] = 0, l<j<m, (2.33) 

as it must be for conserved momenta in quantum mechanics and the Hamiltonian has to 
satisfy the relations 

L r H = , (2.34) 
in order to obtain a Noether symmetry. The vector T is defined by JTTJ 
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Let us now go to the minisuperspace quantum cosmology and to the semi-classical interpre- 
tation of the wave function of the universe. 

By a straightforward canonical quantization procedure, we have 

iij — ► 7tj = —idj , (2.36) 
H^H(q j ,-id qj ). (2.37) 

It is well know that the Hamiltonian constraint gives the WDW equation, so that if > is 
a state of the system (i.e. the wave function of the universe), dynamics is given by 

H\^>=0. (2.38) 

If a Noether symmetry exists, the reduction procedure outlined above can be applied and 
then, from (ggp and ( ggp , we get 



_ 9C 

7Ti = — — = Ixi vr = Z^i , 

dQ 1 1 
dC 

7r a = ^ = <jfA = Sa, (2-39) 
dQ 2 

depending on the number of Noether symmetries. After quantization, we get 

- idx\V > = E x |* > , 

-id 2 \m > = S 2 |^ > , (2.40) 



which are nothing else but translations along the Qi axis singled out by corresponding 
symmetry. Eqs. (|2.40 ) can be immediately integrated and, being Ej real constants, we 



obtain oscillatory behaviours for \ty > in the directions of symmetries, i.e. 

m 

|* >=£ ^ jQi \x(Q l ) > , m<l<n, (2.41) 



=i 



where m is the number of symmetries, / are the directions where symmetries do not exist, n 
is the total dimension of minisuperspace. It is worthwhile to note that the component \x > 
of the wave function could also depend on Sj but it is not possible to state "in general" if 
it is oscillating. 

Viceversa, dynamics given by (|2.38| ) can be reduced by (|2.40| ) if and only if it is possible 



to define constant conjugate momenta as in ( |2.39| ), that is oscillatory behaviours of a subset 
of solutions \ty > exist only if Noether symmetry exists for dynamics. 

The m symmetries give first integrals of motion and then the possibility to select clas- 
sical trajectories. In one and two-dimensional minisuperspaces, the existence of a Noether 
symmetry allows the complete solution of the problem and to get the full semi-classical limit 
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of minisuperspace quantum cosmology. By these arguments, the Halliwell request that an 
oscillatory wave function predict correlations between coordinates and canonical conjugate 
momenta |9j is fully recovered. 

In conclusion, we can set out the following 



Theorem: In the semi-classical limit of quantum cosmology and in the framework of min- 
isuperspace approximation, the reduction procedure of dynamics, due to the existence of 
Noether symmetries, allows to select a subset of the solution of WDW equation where oscil- 
latory behaviours are found. As consequence, correlations between coordinates and canonical 
conjugate momenta emerge so that classical cosmological solutions can be recovered. Vicev- 
ersa, if a subset of the solution of WDW equation has an oscillatory behaviour, due to 
Eq. ( $-4Q ), conserved momenta have to exist and Noether symmetries are present. In other 
words, Noether symmetries select classical universes. 



In what follows, we shall give realizations of such a statement for minisuperspace cosmolog- 
ical models derived from extended gravity theories. 



III. SCALAR TENSOR GRAVITY COSMOLOGIES 



Let us take into account a nonminimally coupled theory of gravity of the form 



A 



J d^x^f^g 



F(<p)R+-g*'<p ll <p v -V(<p) 



(3.1) 



where F(<p) and V(cp) are respectively the coupling and the potential of a scalar field |T3 



We are using, from now on, physical units 8ttG — c — U — 1, so that the standard Einstein 
coupling is recovered for F(cp) = —1/2. 

Let us restrict, for the sake of simplicity, to a FRW cosmology. The Lagrangian in ( |3 . 1| ) 
becomes 



C = 6aa 2 F + 6a z aF — 6kaF + a 



2 ;. 



(3.2) 



in terms of the scale factor a. 

The configuration space of such a Lagrangian is Q = {a, if}, i.e. a two-dimensional 
minisuperspace. A Noether symmetry exists if (|2.15|) holds. In this case, it has to be 



d _ d . d ■ d 
X = a — + p— + a — + P — 

oa 0(f oa oip 



(3.3) 



where a, (3 depend on a, ip. The system of partial differential equation given by ( [2.15|) is 
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F{fp) 



a + 2a 



da 
da 



+ aF'(cp) 



P + a 



dp' 
da 



. , . da d(3 
3a + 12F' {ip — + 2a-^ = 

dip dip 



da dp 
2a + a— + — 

da dip 



apF"(ip) + 

o 

[3aV(<p) + aPV'(<p)}a 2 + Qk[aF(<p) + aPF'( V )} = 



i j i da ^, ^ a 2 dp 

F '( v) + 2-F( v ) + -£ = , 



(3.4) 

(3.5) 

(3.6) 
(3.7) 



The prime indicates the derivative with respect to (p. The number of equations is 4 as it has 
to be, being n = 2. Several solutions exist for this system ]T3}-[r5|. They determine also the 
form of the model since the system (|3.4j) - fl3~77|) gives a,P, F(ip) and V((p). For example, if 
the spatial curvature is k = 0, a solution is 



« = --p(s)/3 a s+ V m(s)_1 , P = A><r ; 



m(s) 



(3.8) 



2p(s) 



where 



D(s) 



F{<p) = D( 8 )<f? , V( V ) = \y 



(2s + 3) 2 . , 3(8 + 1) , s 2s 2 + 6s + 3 
p(s) = , m(s) = 



48(s + l)(s + 2) 



2s + 3 



2s + 3 



and s, A are free parameters. The change of variables ( 2.23 ) gives 

3 



w = CT aV p(s) , z 



where ctq is an integration constant and 



X(s) 



PoX(s 



6s 
2s + 3 



-a ip 



l—m(s) 



Lagrangian ( |3.2|) becomes, for k = 0, 

C = / y(s)w s ^ 3 zw — Xw 

where z is cyclic and 

2s + 3 



70) 



12a 2 (s + 2) (s + 1) ' 



The conjugate momenta are 

dC 

dz 



7T- = = ^f(s)w s/3 W , TT W = ^ = 7(s)w s/3 i . 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



and the Hamiltonian is 



12 



n 



7(s)w s / 3 



+ Aw. 



(3.16) 



The Noether symmetry is given by 



71", 



-*0 • 



Quantizing Eqs. ( |3.15| ), we have 

7r — ► -id z , ir w — ► -id w , 
and then the WDW equation 

[(td z )(td w ) + \w 1+s ^>=o. 

where A = 7(s)A. 

The quantum version of constraint ( ^.17| ) is 

- id z \^ >= Snl^ > , 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



so that dynamics results reduced. A straightforward integration of Eqs. (|3.19| ) and ( |3.20| ) 
gives 



|tf >= \Sl(w) > \x(z) >oc e^° 2 e 



iS z c -i\w 2 + 3 / 3 



(3.21) 



which is an oscillating wave function. In the semi-classical limit, we have two first integrals 
of motion: E (i.e. the equation for n z ) and Ec = 0, i.e. the Hamiltonian ( 3.16| ) which 
becomes the equation for n w . Classical trajectories in the configuration space Q = {w , z} 
are immediately recovered 



w(t) = [ht + hf^ , 
z(t) = [k 1 t + k 2 ]^^ s+ ^ + z , 

then, going back to Q = {a, <p}, we get the classical cosmological behaviour 

a(t) = a (t-t ) 1 ^, 



where 



1(8) 



2s 2 + 9s + 6 
s(s + 3) 



q(s) 



2s + 3 



(3.22) 
(3.23) 



(3.24) 
(3.25) 



(3.26) 



Depending on the value of s, we get Friedman, power-law, or pole-like behaviours. 

If we take into account generic Bianchi models, the configuration space is Q = 
{ai, 02, az, (p} and more than one symmetry can exist as it is shown in |12| . The con- 
siderations on the oscillatory regime of the wave function of the universe and the recovering 
of classical behaviours are exactly the same. 
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IV. FOURTH-ORDER GRAVITY COSMOLOGIES 



Similar arguments work for higher-order gravity cosmology. In particular, let us consider 
fourth-order gravity given by the action 



A = J d A x^=gf(R) 



(4.1) 



where f(R) is a generic function of scalar curvature. If f(R) = R+2A, the standard second- 
order gravity is recovered. We are discarding matter contributions. Reducing the action to 
a point-like, FRW one, we have to write 



A = / dt£(a, a; R, R) 



(4.2) 



where dot means derivative with respect to the cosmic time. The scale factor a and the Ricci 
scalar R are the canonical variables. This position could seem arbitrary since R depends 
on a, a, a, but it is generally used in canonical quantization [IB-TS]. The definition of R in 
terms of a, a, a introduces a constraint which eliminates second and higher order derivatives 
in action (|4.2j ), and yields to a system of second order differential equations in {a, R}. Action 
L2I) can be written as 



A = 2ir 2 J dtla 3 f{R) - A 



fa a k 
R + 61- + — + — 



(4.3) 



where the Lagrange multiplier A is derived by varying with respect to R. It is 

A = a 3 f'(R) . (4.4) 

Here prime means derivative with respect to R. To recover a more strict analogy with 
previous scalar-tensor models, let us introduce the auxiliary field 



p = f(R) 



(4.5) 



so that the Lagrangian in (43) becomes 

C = 6aa 2 p + 6a 2 dp — 6kap — a 3 W(p) 



(4.6) 

which is of the same form of (|3.2| ) a part the kinetic term. This is an Helmhotz-like La- 
grangian |BJ and a,p are independent fields. The potential W(p) is defined as 



W(p) = h{p)p — r(p) 



where 



r(p) = J f(R)dR = JpdR = f(R) , h(p) = R. 



(4.7) 



(4.8) 



such that h = (/') 1 is the inverse function of /'. The configuration space is now Q = {a,p} 
and p has the same role of the above ip. Condition ( J2.15 ) is now realized by the vector field 
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d d d • d 

X = a(a,p) — + P(a,p)— + a— + (3— 
da dp da op 



(4.9) 



and explicitly it gives the system 



P 



a + 2a 



da 
da 



p + a 



13 + a 



dp 
da 



9 da 
a 2 — = 0, 
dp 

da da d(3 
2a + a— + 2p— + a^f = . 
da dp dp 



6k[ap + Pa] + a 2 [3aW + af3 



dW 
dp 



0. 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



The solution of this system, i.e. the existence of a Noether symmetry, gives a, (3 and W(p). 
It is satisfied for 



a = a(a) , (3(a,p) = /3 a s p, 
where s is a parameter and (3q is an integration constant. In particular, 

s = -^ a (a) = -^a, f3(p) = p oPl W(p) = W oP , k = 



(4.14) 



(4.15) 



a (a) 



p(a,p) = (3 P-, W(p) = W lP 3 , V k. 



(4.16) 



where Wq and W\ are constants. As above, the new set of variables Qi = Q\q % ) adapted 
to the foliation induced by X are given by the system ( |2.23 ). Let us discuss separately the 
solutions ( CT ) and (^T6|) . 



A. The case s = 



The induced change of variables Q = {a,p} — >■ Q = {w, z} can be 

w(a,p)=a 3 p, z(p)=\np. 

Lagrangian ( f4.6| ) becomes 



if 



£(w, w, z) = zw — 2wz H 3Wow . 

w 

and, obviously, z is the cyclic variable. The conjugate momenta are 

7T Z = T7T = W ~ AZ = S , 



(4.17) 



(4.18) 



(4.19) 
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and the Hamiltonian is 



7T 

k w ti z + 2wnl + 6W w . 

w 



By canonical quantization, reduced dynamics is given by 



d' z z - 2w 2 dl - wd w d z + 6W w' 



Itf >= 



(4.21) 



(4.22) 
(4.23) 



However, we have done simple factor ordering considerations in the WDW equation ( |4.22[ ). 
Immediately, the wave function has an oscillatory factor, being 



The function \x > satisfies the Bessel differential equation 



w 2 dl + i^wd w +(^-3W Q w* 



X(w) =0 



whose solutions are linear combinations of Bessel functions Z v {w) 

x (w) = w l / 2 ^l A Z v (\w) , 

where 

1 



v = ±- 



4 - 9Sg - i4S c 



A = ±9i 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



4 v u y 2 

The oscillatory regime for this component depends on the reality of z/ and A. The wave 
function of the universe, from Noether symmetry ( 4.15|) is then 

V(z,w) ~ e iEo[2 - (1/4)lnwl w 1/2 Z u (Xw) . (4.28) 



For large w, the Bessel functions have an exponential behaviour pOfl , so that the wave 
function ( 4.28| ) can be written as 

o i[S 2-(S /4) \nw±Xw] 



(4.29) 



By identifying the exponential factor of Q4.2S ) with So, we can recover the conserved mo- 
menta 7r z , tc w and select classical trajectories. Going back to the old variables, we get the 
cosmological solutions 



a(t) = a e^ exp \-^- e^ 2 ^ 



p(t)=p Q e^ t exp{^ ie -( 2A /3)*} 



(4.30) 



(4.31) 



where do, Po and z\ are integration constants. It is clear that A plays the role of a cosmological 
constant and inflationary behaviour is asymptotically recovered. 
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B. The case s = -2 
The new variables adapted to the foliation for the solution (|4.16|) are now 



w(a,p)=ap, z(a) = a 2 . (4.32) 
and Lagrangian ( |4.6|) assumes the form 

C(w, w, z) = 3zw — 6kw - WiW 3 , (4.33) 
The conjugate momenta are 

^ = — = 3^ = ^, (4.34) 
<9£ 

k w = — = 3z. (4.35) 
ow 

The Hamiltonian is given by 

H(w, tt w , 7r 2 ) = - n z n w + 6kw + Wxw 3 . (4.36) 
Going over the same steps as above, the wave function of the universe is given by 

and the classical cosmological solutions are 

a(t) = ±fi®, p{t) = ± C _1±^B1 , (4.38) 

where 

h(t) = (m.^ «. + (p^j f + ( Kl + e + W1 (6A + t + , 2 . (4.39) 

Wi, Z\ and z 2 are integration constants. Immediately we see that, for large t 

a(t) ~ t 2 , p(t) ~ - . (4.40) 
which is a power-law inflationary behaviour. 
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V. HIGHER THAN FOURTH ORDER GRAVITY COSMOLOGIES 



Minisuperspaces which are suitable for the above analysis can be found for higher than 
fourth-order theories of gravity as 

A = J <Pxy/=gf(R,DR). (5.1) 

In this case, the configuration space is Q = {a, R, OR} considering OR as an independent 



degree of freedom |p!7|, ^8| , pl]| . The FRW point-like Lagrangian is formally 

C = C(a,a,R, R,OR,(OR)) 

and the constraints 



R=-6 



a / a\ 2 k 
-+(-) +- 



(5.2) 



(5.3) 



OR = R + 3-R 

a 



(5.4) 



holds. Using the above Lagrange multiplier approach, we get the Helmholtz point-like 
Lagrangian 



where 



£ = Qaa p + 6a dp — 6kap — a hq — a W(p, q) 



df 



df 



P ~ OR ' q dOR 



(5.5) 



(5.6) 



and 



W(p, q) = h(p)p + g(q)q - f 



(5.7) 



h(p) = R, g(q) = OR, f = f(R,OR). (5.8) 
Now the minisuperspace is three-dimensional but, again, the Noether symmetries can be 



recovered. Cases of physical interest [21 



are 



f(R, OR) = F R + F 1 R 2 + F 2 ROR . 



(5.9) 



f(R,OR) = F R+F 1 VROR. 



(5.10) 



discussed in details in [p2]| . Also here the existence of the symmetry selects the form of the 
model and allows to reduce the dynamics. Once it is identified, we can perform the change 
of variables induced by foliation using Eqs. (|2.23|) , if a symmetry is present, or Eqs. ( p.25|) , 
is two symmetries are present. In both cases, 
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Q = {a, R, OR} — ► Q = {z, u, w} 



(5.11) 



where one or two variables are cyclic in Lagrangian Q5.5| ). Taking into account, for example, 
the case (|5.10| ), we get 



C = 3[ww 2 — kw] — Fi 



3ww 2 u + 3w 2 wu + 



w 3 zii 



2u 2 



3kwu 



(5.12) 



where we assume F = —1/2, the standard Einstein coupling, z is the cyclic variable and 

(5.13) 



z = R, u 



R 



w = a . 



The conserved quantity is 



w 3 u 



2u 



2 ' 



(5.14) 



Using the canonical procedure of quantization and deriving the WDW equation from ( p. 12 ), 
the wave function of the universe is 



>^ e iEoZ \ X {u) > \G{w) > 



(5.15) 



where xi u ) an d Q{w) are combinations of Bessel functions. The oscillatory subset of the 
solution is evident. In the semi-classical limit, using the conserved momentum fl5.14 ), we 
obtain the cosmological behaviours 



a(t) = dot , a(t) = aot 1 ^ 2 , a(t) = a Q e kot . 
depending on the choice of boundary conditions. 



(5.16) 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we discussed the connection of Noether symmetries for minisuperspace 
cosmological models to the recovering of classical solutions. If the wave function of the 
universe is related to the probability to get a given classical cosmology, the existence of such 
symmetries tell us when the WDW wave function of the universe has oscillatory behaviours 
connected to the recovering of correlations between coordinates and conjugate canonical 
momenta . I this sense, the Hartle criterion to get correlations capable of selecting classical 
universes works. 

Some remarks are necessary at this point. First of all, we have to stress that the wave 
function is only related to the probability to get a certain behaviour but it is not the prob- 
ability amplitude since, till now, quantum cosmology is not a unitary theory. Furthermore, 
the Hartle criterion works in the context of an Everett-type interpretation of quantum cos- 
mology |23J23] which assumes the ideas that the universe branches into a large number of 



copies of itself whenever a measurement is made. This point of view is called Many Worlds 
interpretation of quantum cosmology. Such an interpretation is just one way of thinking 
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and gives a formulation of quantum mechanics designed to deal with correlations internal 
to individual, isolated systems. The Hartle criterion gives an operative interpretation of 
such correlations. In particular, if the wave function is strongly peaked in some region of 
configuration space, we predict that we will observe the correlations which characterize that 
region. On the other hand, if the wave function is smooth in some region, we predict that 
correlations which characterize that region are precluded to the observations. 

If the wave function is neither peaked nor smooth, no predictions are possible from 
observations. In other words, we can read the correlations of some region of minisuperspace 
as causal connections. However, the validity of minisuperspace approximation is often not 
completely accepted and it is still matter of debate . 

As we said above, the analogy with standard quantum mechanics is straightforward. By 
considering the case in which the individual system consists of a large number of identical 
subsystems, one can derive from the above interpretation, the usual probabilistic interpre- 
tation of quantum mechanics for the subsystems |T|Q]. 

What we proposed in this paper is a criterion by which the Hartle point of view can be 
recovered without arbitrariness. If a Noether symmetry (or more than one) is present for 
a given minisuperspace model, then oscillatory subsets of the wave function of the universe 
are found. Viceversa, oscillatory parts of the wave function can be always connected to 
conserved momenta and then to Noether symmetries. 

^From a general point of view, this is the same philosophy of many branches of physics: 
Finding symmetries allows to solve dynamics, gives the main features of systems and simplify 
the interpretation of results. 

However the above scheme should be enlarged to more general classes of minisuperspaces 
in order to seek for its application to the full field theory (i.e. to the infinite-dimensional 
superspace). Only in this sense, one could claim for the validity of the approach to the full 
semiclassical limit of quantum cosmology. 
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